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A TERRIBLE EXPANSION OF THE DETERMINANT 


ERIK INSKO, KATIE JOHNSON, AND SHAUN SULLIVAN 


Abstract. From a transfer formula in multivariate finite operator calculus, comes an expansion 
for the determinant similar to Ryser’s formula for the permanent. Although this one contains 
many more terms than the usual determinant formula. To prove it, we consider the poset of 
ordered partitions, properties of the permutahedron, and some good old fashioned combinatorial 
techniques. 


1. Introduction 


One of the foundational concepts of linear algebra is the determinant. At the most basic level, 
this matrix parameter is celebrated for its intricate ties to the set of eigenvalues and as a similarity 
invariant. However, the determinant still surprises us as the solution to a varying array of problems. 

In addition to solving systems of linear equations and performing a change of variables in calculus, 
the determinant can help us count! Benjamin and Cameron [1] recently showed the determinant 
will calculate the number of nonintersecting n-paths in certain nonpermutable digraphs, where an 
n-path is a set of n paths from n distinct source vertices to n distinct sink vertices. In fact, the 
permanent will count the number of all n-paths. 

The determinant of a matrix can be found recursively, as an alternating sum of minors. Often 
the determinant of an n x n matrix A is defined compactly using the Leibniz formula, precisely 

n 

det(A) = ^ sgn(cT) 

O'&Sn *=1 

Similarly, the permanent can be defined as a sum over subsets of [n] := {1, 2,..., n} using Ryser’s 
formula [3] 

n 

perm(A) = X 11 X 

Sc[n] i=lj£S 

In this paper, we prove a much messier expansion of the determinant by instead indexing our 
terms using the set of ordered partitions of [n]. Aptly, we call this the terrible expansion of the 
determinant. This expansion is analogous to Ryser’s formula for the permanent. Section [S] explains 
the origins of this expansion as it relates to multivariate finite operator calculus, a branch of 
mathematics that has proven useful in enumerating ballot (generalized Dyck) paths containing 
certain patterns 017]. 

Before stating the formula for our expansion of the determinant, we introduce it with two fun¬ 
damental examples. When n = 2 we see that Equation ([T]) provides the following expansion for the 
determinant. 


( 1 ) 


ail cli 2 
021 0,22 


aii(ai 2 + 022) + 022(011 -I- 021) — (on -|- a 2 i)(oi 2 + 022) 
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Likewise when n = 3 we get the following expansion: 

1^1 = ail (ai2 + 022) (ai3 + 023 + 033) + 011(013 + 033) (ai2 + 022 + 032) 

+022(011 + a2i)(ai3 + 023 + 033) + 022(023 + a33)(oii + 021 + 031) 

+033(011 + a3i)(ai2 + 022 + 032) + 033(022 + a32)(oii + 021 + 031) 

— 011(012 + O22 + 032)(ai3 + 023 + O33) — 022(011 + 021 + a3l)(oi3 + O23 + O33) 

— 033(011 + O21 + 03i)(oi2 + O22 + O32) — (oil + a2l)(ai2 + 022)(oi3 + 023 + O33) 

— (oil + 03l)(0i3 + a33)(oi2 + 022 + O32) — (022 + a 32 ){a 23 + 033)(aii + O21 + 031) 
+ (oii + O21 + a3i)(oi2 + 022 + 032)(ai3 + O23 + 033). 


Our main theorem gives a general description of the terrible expansion of the determinant. 
Theorem 1. Let A = (ojj)^^^. The following formula is an expansion for the determinant of A: 

(2) det(yi) = (-i)”-i®i n n E “«• 

where the outer summation runs over all ordered partitions B = (/3i,/32, ■ ■ ■, Pr) of the set [n] and 
the inner summation runs over all integers i in the union of first k parts fj'f, = U^=i l^j of the 
partition B. 

The following example serves to clarify the notation in Theorem [TJ 


Example 2. If 

i 3 = (/ 3 i,/ 32 ,/? 3 ) = ({ 2 },{l, 3 },{ 4 , 5 }), 

then 

B' = (/3(,/3',/3') = ({2}, {1,2,3},{1,2,3,4,5}) , 
and the corresponding expression in Equation ([2]) for the ordered partition B is 

022(011 + O21 + a3i)(oi3 + O23 + a33)(oi4 + 024 + 034 + O44 + 054)(oi5 + O25 + O35 + O45 + O55). 
notice that the first index runs through B', while the second index runs through the partition B. 

The rest of this paper proceeds as follows; In Section [2l we analyze the functions / : [n] —>■ [n] 
indexing the terms in the terrible expansion. After setting notation and proving a few fundamental 
lemmas, we end that section with Corollary 1101 which proves that when / : [n] —>■ [n] is a bijective 
function, or a permutation, then the coefficients Cf = sgn{f). In other words, they are precisely 
the nonzero coefficients appearing in the determinant. In Section O we study the poset of ordered 
partitions and identify the importance of singleton partitions so that we can formulate our problem 
in more geometric terms as Euler characteristics of convex polytopes relating to the permutahedron. 
In Section m we prove that cj = 0 for all non-bijective functions / : [n] — )■ [n] by analyzing Euler 
characteristics of subsets of the permutahedron. This proves that the terrible expansion does 
indeed give a formula for the determinant. In Section [ 5 l we give an extremely brief introduction to 
multivariate finite operator calculus, and state a more general open conjecture that motivated this 
paper. 


2 . Flattening Functions 

Upon expanding the expression in Theorem [H many terms will cancel. In this section, we set up 
the groundwork to keep track of each of the terms and show how they cancel. Now consider any 
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function / : [n] —>■ [n] and define the monomial 

n 

Expanding the terms in Equation ([2|) resnlts in a sum of the form 

(3) E n n 

Bh[n] 0k&B f 

where each term CfUf corresponds to a function from the set [n] to itself. Eor instance, in Equation 
© the four functions /* : [ 2 ] [ 2 ] are 

/i(l) = 1 , /i( 2 ) = 1 ; / 2 ( 1 ) = 1 , / 2 ( 2 ) = 2 ; / 3 ( 1 ) = 2 , ^(2) = 2 ; and Ml) = 2 , M‘^) = 1 - 
The terms corresponding to each function are labeled below 

aii(ai2 + 022) + 022(011 + 021) —(oil + 02i)(ai2 + 022) 

= (011012 — 011012) + (011O22 + 011022 — 011022)+ 

(021022 — 021022) + (—012021) 

=0 + 011022 + 0 — O12O21 

=C/lO/i + C/ 2 O /2 + C /3 0/3 + C /4 0 / 4 , 

where first we expand the terms and then we simplify. Hence we see that c/^ = 0, c /2 = 1, c /3 = 0, 
and c /4 = — 1 . 

Thus, the goal of this paper is to combinatorially identify the coefficients c/ for each such function, 
and show that they agree with the coefficients of o/ in the determinant. To do so, we must identify 
the set Sf = {B h [n] : af appears as a summand of the product indexed by B}. The following 
definition and lemma describe criteria for when an ordered partition B b [n] appears in S'/. 

Definition 3. Let B = (/?i, /32, • • •, /^r) be an ordered partition of [n]. For each i G [n], define 
/3(i) = A: if i G /3fc. We say i precequals j in B, denoted i + j, if i appears in an earlier part or 
the same part as j in the ordered partition B, i.e., 

m < m- 

Lemma 4. Let B h [n] he an ordered partition of [n]. The term a/ appears in the product 
n n E iff fU) :<j in B for all 1 < j < n. 

Proof. Let B h [n] be an ordered partition of [n]. Since j3'^ = U/=i/3/, we see that a/ appears as a 
term in the product 0 0 E ^ij precisely when i = /(/)+ j in the ordered partition B for all 

1 < j <n. □ 

Now that Lemma m has established that an ordered partition B appears in the set S'/ if and 
only if f{j) precequals j in B for all integers 1 < j < n, we are ready to analyze which functions 
/ : [n] —>■ [n] correspond to nonzero terms in the terrible expansion. To do so, we start by 
introducing some notation regarding the structure of such functions. 

Definition 5. Let f : [n] ^ [n] be any function. We say that the function / is acyclic if f^{i) = i 
implies A: = 1 for all i G [n]. Otherwise, we say / contains a cycle . 

Henceforth we will prefer to work with acyclic functions. The next definition, describes how each 
function / : [n] —> [n] which contains a cycle can be simplified or flattened to an acyclic function / 
on a different set of elements. 
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Definition 6. Let / : [n] — >■ [n] be any function, Cf = {cti, (T2, ..., (t^} represent the cycles of 
/; Nf C [n] be the elements not in a cycle, and Df = Cf Li Nf. Define the flattened function 
f : Df ^ Df as follows 


/(*) 


f{i) if f{i) is not in a cycle of / 

< aj if f{i) belongs to the cycle aj ■ 
i if i Cf 


Intuitively, / acts just like /, but shrinks each cycle of / to a fixed point, and thus it is an acyclic 
function. The following example illustrates Definitions [5] and [6l 


Example 7. Let / : [6] ^ [6] be defined by /(I) = l,/(2) = 3,/(3) = 2, /(4) = 3, /( 5 ) = 6 , 
and /(6) = 5 . Then Cf = {fTi,cT 2 } where ui = (23) and C 2 = ( 56 ). The sets Nf = {1,4} and 
Df = Nf L) Cf = {(Ti,(T 2 , 1,4}. The function / ; Df Df is defined by /(ui) = ui, /(cr 2 ) = (T 2 , 
/(I) = 1, and 7(4) = ui. 

The following lemma shows that we can reduce the problem of calculating the coefficients Cf in 
the terrible expansion, to that of calculating the coefficients Cf corresponding to acyclic functions. 

Lemma 8. If f : [n] ^ [n] is any function, then the coefficients Cf and cj are related by the 
equation Cf = (—l)”“l'^/lcj. 

Proof. Let Sf := {B h [n] : Of appears as a summand of the product indexed by B}, and similarly 
let Sj := {A h Df : aj appears as a summand of the product indexed by A}. 

If B is an ordered partition for which Of appears as a summand, then Lemma H] implies that for 
each cycle 

of /, the following precedence relation must hold 

i^f\i)^f'^-\i)^---^f{if^i 

in the ordered partition B. Therefore Cj must be contained in the same part /3 of the ordered 
partition B. From here, it is clear to see that Sf is equivalent to the set of ordered partitions of 
Df, and so each term in Equation ([3]) has the form 

c,a, = Y, 

B&Sf 

A&Df 

and the lemma follows. □ 


The number of ordered partitions of [n] with k parts is well known to be k\S{n, k), where S{n, k) 
are the Stirling numbers of the second kind. We obtain an important corollary to Lemma [5] from 
the following well-known result about Stirling numbers of the second kind. 

Lemma 9. The following identity holds for ordered partitions: 

n 

^{-ir-’^k\S{n,k) = l. 

k=0 
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Proof. The result follows immediately upon setting x = —1 in the following identity on Stirling 
numbers of the second kind[ 6 l p. 35]: 

n 

^^S{n,k){x)k = x^. □ 

k=0 

Corollary 10. If f : [n] — 7 > [n] is bijective, i.e., / = vr for some vr G then Cf = sgn( 7 r). 

Proof. Since / is bijective, it consists only of cycles. Thus Dj = Cf = {ai,... ,ar} and \Df \ = r. 
By Lemma [ 8 ] and Lemma [9] we see 

Aev 

B\-[r] 

n 

k=0 

= (-l)^-L 

We leave it to the reader to verify that sgn( 7 r) = (—1)’^“'’. □ 

With Corollary [TOl we have that the terrible expansion contains every term of the determinant 
with the correct coefficient. It remains to show that whenever / is not bijective, Cf = 0 . By 
Lemma [8l it suffices to show this for acyclic functions. 

3 . The Poset of Ordered Partitions 

We next consider the poset of ordered partitions in order to show that the set Sf has a nice 
structure when / is acyclic. This will allow us to eventually switch to a more geometric viewpoint. 

Let Vn denote the poset of ordered partitions of the set [n]. In Figure [U we see the poset P 3 
of ordered partitions on 3 elements. At the top of the poset we have the n! ordered partitions 
consisting of singletons. These partitions correspond bijectively with the elements of Gn- Because 
of their importance later, we will refer to them as singleton partitions. 



Figure 1. Poset of Ordered Partitions on 3 elements 

Directly below a given ordered partition B in Vn are ordered partitions formed by taking the 
union of two consecutive parts in B. For example, directly below the singleton partition 3/1/2/A 
are the ordered partitions 13/2/4, 3/12/4, and 3/1/24. All ordered partitions under a singleton 
partition creates an (n — l)-cube. An example is given in Figure [2j 
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An acyclic function can be viewed as a rooted forest, where the fixed points are the roots. An 
example is given in Figure [3l Given an acyclic function /, if a path exists from p to q, with p closer 
to the root than q, then f^{q) = p for some k. Thus, p ^ q, and so we say / has the rule p ^ q. In 
this way, each function / stipulates a set of rules 

Rf = {p^q \ f^{q) =pioT q,ke [n]| . 

The following lemma and corollary will show that the set of ordered partitions Sf has a nice 
structure in Vn- 


1/2/3/4 



1/2/34 


1/234 


Figure 2. The Cube Under the Singleton Partition 1/2/3/4 



Figure 3. Acyclic Function Represented by a Rooted Forest 


Lemma 11. If a singleton partition A in Vn satisfies the rules Rf given by an acyclic function f 
then every ordered partition B < A in Vn also satisfies the rules Rf. If an ordered partition B in 
Vn satisfies the rules Rf then there exists at least one singleton partition A> B inVn that satisfies 
the precedence rules Rf. 

Proof. The first statement is obvious. If A is a singleton partition with the rule p P q in A, and 
R < A is an ordered partition in Vn, then the parts of B are unions of consecutive parts of A. 
Hence p ^ g in R as well. 

The second statement is slightly less obvious. Let / be an acyclic function, and suppose that 
R in Vn is a nonsingleton ordered partition that satisfies the rules Rf. We must show there is a 
singleton partition A > B above it in Vn that also satisfies Rf. Consider a part fi of B that is not 
a singleton. If no pair of elements in fi has a rule associated with it, then the elements of j3 can be 
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ordered arbitrarily. Otherwise, there are elements of (3 that have rules imposed on them. Consider 
the elements in the intersection of {3 and a rooted tree associated with /. We order those elements 
by their distance from the root. (Those elements having the same distance from the root can be 
put in any order with respect to each other.) We do this for every rooted tree associated with / 
to impose an order on all of /?. Doing the same to each part will result in a singleton partition A 
above B satisfying the rules of /. □ 

Corollary 12. The set of ordered partitions of [n] satisfying the rules of an acyclic function is a 
union of (n — 1) -cubes in Vn- 

Proof. By the above lemma, we can account for all the ordered partitions by only considering the 
singleton partitions satisfying the acyclic function, and all the ordered partitions below them. The 
result follows since the ordered partitions below a singleton partition form an (n — l)-cube. □ 

Lemma [H] and Corollary [12] tell us that once we know which singleton partitions appear in S'/, 
then we know that S/ is precisely those singletons and all the ordered partitions under them in 
Vn- Because of their importance we will start to label the singleton partitions without slashes, e.g. 
1/2/3/4 —>■ 1234, unless we need to distinguish them from the ordered partition with one part. In 
the next section, we turn our attention to a geometric object isomorphic to Vn, the permutahedron. 

4 . The Permutahedron 

In this section we show how the alternating sums giving c/ when / is acyclic are related to 
the Euler characteristic of the permutahedron and use this correspondence to show that c/ = 0 . 
It is well-known that the poset of the ordered partitions is isomorphic to the face lattice of the 
permutahedron [5] Fact 4.1]. Specifically, each vertex on the permutahedron represents a singleton 
partition, the edges incident to a vertex represent the ordered partitions just below that singleton 
partition in the poset, the faces adjacent to those edges represent the ordered partitions just below 
again, and so on, until the permutahedron itself represents the ordered partition with one part at 
the bottom of the poset. Note that two vertices are adjacent if one can be obtained by a single 
swap of consecutive elements. For example, 315624 is adjacent to 351624. Figure ID shows the 
transformation from the poset on 3 element to the permutahedron on 3 elements, which in this case 
is a hexagon. Figure [5] shows the poset on 4 elements as the ordinary permutahedron (truncated 
octahedron). 


1/2/3 



1/2/3 



Figure 4. Poset of Ordered Partitions on 3 elements: Top View (left) and as a 
Permutahedron (right) 


7 









4321 



Figure 5. Poset of Ordered Partitions on 4 elements as a Permutahedron 

In Figure [SJ we only label the singleton partitions at the vertices, but the labeling of the other 
ordered partitions would be similar to Figure [H 

The permutahedron is often defined as the convex hull of the points 

Pa = (ct(1),(t(2), .. .,a{n)) 

for every a G 6„. It is a convex polytope, and in particular, it is contractible to a point. Thus, 
the permutahedron has Euler characteristic 1 [5]. Since the ordered partitions are in bijection with 
the faces of the permutahedron, the alternating sum of the ordered partitions is precisely the Euler 
characteristic of the permutahedron, and this gives us a second proof of Lemma [9l 

We adopt a slightly different convention, relabeling the vertices of the n„ to We will 

denote this relabeled permutahedron by 11^ . Eigure [6] shows IIs in the relabeled 11^, and the 
correspondence between xi < X 2 and 1^2. Figure [6] also shows the fact that n„, and thus 11^, is 

an (re —l)-dimensional object, since all the points lie in the hyperplane Xi + X 2 + - • ■+Xn = 

In general Xi < Xj in n„ corresponds to i ^ j in 

Lemma 13. The singleton partitions in satisfying the precedence rule i ^ j are contained in 
the corresponding half-space Xi < Xj in Un- 

Proof. A permutation satisfies ^ precisely when a{i) < cr{j). Hence we see that if 
the permutation (or singleton partition) satisfies the precedence relation i ^ j then the vertex 
((t(1), (t( 2), ... , (T(re)) is in the half-space Xi < Xj and vice versa. □ 

We are now ready to prove the main result of this section. 

Proposition 14. If f : [re] ^ [re] is not bijective, then the coefficient Cf = 0 is zero. 

Proof. Let / : [re] —[re] be an acyclic function. Let Rj = {/(j) = i ^ j} denote the set of 
precedence rules determined by /. To each precedence rule /(j) = i :< j in Rf we can assign a 
half-space Hij := {(xi, X2 ,... ,Xn) G K"" : = Xi < Xj}, and the intersection of these half-spaces 

with the permutahedron H]^ defines a convex polytope, which we will denote by Hy. 

The faces of Hj fall into two disjoint subsets: those faces that correspond to the ordered partitions 
in Sf and those that do not. Let T f denote the faces of Hj which correspond to elements of S'/ 
and let A/ denote the faces of H/ that do not. The set of faces A/ are precisely the faces of H/ 
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Figure 6 . The permutahedra Ha (left) and Ilg (right) 


which lie entirely on the boundary of at least one half-space Xi = Xj because they resulted from 
intersecting with one of the half-spaces Hij. Thus each face of Aj is a convex polytope, and 
we see that Ay is a union convex polytopes. Each of these convex polytopes contains the point 
(xi, X 2 ,..., Xn) where xi = X 2 = • • • = x^, so Ay is a contractible space. Hence Ay has Euler 
characteristic xi^f) = 1 - 

Since Hy is a convex poly tope and ny=ryLlAywe see that 

1 = x(ny) = x(ry) x(Ay) = 0 -h 1. 

It follows that x(ry) = 0. Since the summands of the alternating sum (_i)’^“l^l correspond 

BeSf 

with the faces of Ey we see that cy = (_i)’^-l^l = ^(py) = 0 as desired. □ 

BeSf 


We end this section with an example demonstrating the proof of Proposition [TTl 


Example 15. Eigure[7]shows Hy with T y bold and Ay shaded for the terms oy = 011012013 , which 
has the rules 1^2 and 1^3 and oy = 011012033034 , which has the rules 1^2 and 3^4. 


5 . Multivariate Einite Operator Calculus 

This terrible expansion of the determinant came from a conjecture about a transfer formula in 
multivariate finite operator calculus (MEOC). In this section we give a very brief overview of the 
objects of study in MEOC and the conjecture that gives this expansion. The interested reader is 
encouraged to read [ 8 ] for a more comprehensive description of this subject matter. 

Let k he a held. Let denote the standard ^-dimensional basis of . The main objects 

of study in MEOC are polynomials p G A;[xi, ..., x^] and operators T € A:[[i4i, ..., L*r]], where Di 
is the partial derivative with respect to Xj. A sequence of polynomials 6 n(x) = bnx,...,ni{xi, ...^xfi) 
is called a Sheffer sequence if there is a set of operators B = (Hi,... ,Bi) with Bi = DiPi where 
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Figure 7. o/ = 011012013 (left) and aj = 011012033034 (right) 


each Bi : ^ ^n-ei and each Pi an invertible operator. Such a set of operators is called a delta 

£-tuple. The power series for an operator T is written as 


(4) r = ^OnD"= ^ o„,••• 

n>0 

These are all standard notations in any multivariate theory. However, the following notation is 

not completely standard in MFOC. Given a subset A C [P\ we define Xa '■= 0 

i&A 

The following is the Transfer Theorem from MFOC, and is essentially Theorem 1.3.6 in [8] with 
different notation. 

Theorem 16 (Transfer Theorem). Let J denote the usual Jacobian matrix of a collection of 
polynomials. Suppose B = (Hi, H 2 ,..., Bf) is a delta i-tuple where Bi = DiP~^, then 


bJx ) = P 


n+l 


J{Bi,B2, ... ,Bi) — 

n \ 


is the basic sequence for B written in terms of ■ 


dBi 


Within the Jacobian, we have Pincherle derivatives —— = BiOj — OjBi, where 04 : p —>■ Xjp is the 

oDj 

jih. umbral shift operator that does not commute with the delta operators. Thus, there are many 
ways this transfer formula could be expanded. The following conjecture (based on the examples 
provided below) gives one such way. 


Conjecture 17. The basic sequence from the Transfer Theorem can also be calculated as 


= E (-1)'-'®' (Ol>Pf)B 

Bt\ei 

where B runs through all ordered partitions of [^] and fd runs through the partitions of B. 
Because of our abuse of some notation, we give some examples. 


6m,n(R, I^) = {uPf^vP^ + vP^uPf^ - UvPf^P^ 


ym—lyn—1 
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= (aRbScT + aRcTbS + bSaRcT + bScTaR + cTaRhS + cTbSaR 


—aRbcST — bSacRT — cTabRS — abRScT — acRTbS — bcSTaR 


+abcRST) ■ 


B = ({2}, {1,3}) ^ {Oi3P0)b = bSacRT 

The terrible expansion of the determinant comes from setting each Bi = D®* = 
or in other words, it is one term of the power series in Equation (j4]). 

6 . Conclusions and Open Questions 

We end with few open questions stemming from our work. 

Ql. Now that Theorem [1] shows that Conjecture [T7] is true for one term of an operator’s power 
series, can Conjecture [T3 be proven by a linearity argument? 

Q2. Can the proof of Ryser’s formula given by Horn and Johnson [3] be modified to give another 
proof of Theorem [TI’ 

Q3. Can our proof of Theorem [T] be modified to prove Ryser’s formula by using the topologi¬ 
cal/combinatorial properties of the cube instead of the permutahedon? 
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